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Heavy Ion Collisions (HIC) represent a unique tool to probe the in-medium nuclear 
interaction in regions away from saturation. In this work we present a selection of reac- 
tion observables in dissipative collisions particularly sensitive to the isovector part of the 
interaction, i.e. to the symmetry term of the nuclear Equation of State (EoS). At low en- 
ergies the behavior of the symmetry energy around saturation influences dissipation and 
fragment production mechanisms. We will first discuss the recently observed Dynami- 
cal Dipole Radiation, due to a collective neutron-proton oscillation during the charge 
equilibration in fusion and deep-inelastic collisions. Important Iso — EOS are stressed. 
Reactions induced by unstable ^^-^ Sn beams appear to be very promising tools to test 
the sub-saturation Isovector EoS. New Isospin sensitive observables are also presented 
for deep-inelastic, fragmentation collisions and Isospin equilibration measurements (Im- 
balance Ratios). 

The high density symmetry term can be derived from isospin effects on heavy ion 
reactions at relativistic energies (few AGeV range), that can even allow a "direct" study 
of the covariant structure of the isovector interaction in the hadron medium. Rather 
sensitive observables are proposed from collective flows and from pion/kaon produc- 
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tion. The possibility of the transition to a mixed hadron-quark phase, at high baryon 
and isospin density, is finally suggested. Some signatures could come from an expected 
"neutron trapping" effect. The importance of studying violent collisions with radioactive 
beams from low to relativistic energies is finally stressed. 

Keywords: Charge Equilibration, Symmetry Energy, Isospin Transport, Isospin Flows, 
Particle Production, Deconfinement 



1. Introduction 

The symmetry energy Egym appears in the energy density e(p, ^3) = e{p) + 
pEsy7n{p3/ p)^ + 0{p^/pY + •■, expressed in terms of total {p = Pp + Pn) and 
isospin (p3 = Pp — Pn) densities. The symmetry term gets a kinetic contribution 
directly from basic Pauli correlations and a potentialpart from the highly contro- 
versial isospin dependence of the effective interactions^. Both at sub-saturation and 
supra-saturation densities,_predictions based of the existing many-body techniques 
diverge rather widely, see We remind that the knowledge of the EoS of asym- 
metric matter is very important at low densities (neutron skins, nuclear structure 
at the drip lines, neutron distillation in fragmentation, neutron star formation and 
crust..) as well as at high densities (transition to a deconfined phase, neutron star 
mass/radius, cooling, hybrid structure, formation of black holes...). We take ad- 
vantage of new opportunities in theory (development of rather reliable microscopic 
transport codes for HIC) and in experiments (availability of very asymmetric ra- 
dioactive beams, improved possibility of measuring event-by-event correlations) to 
present results that are severely constraining the existing effective interaction mod- 
els. We will discuss dissipative collisions in a wide range of energies, from just above 
the Coulomb barrier up to a few AGeV. Low to Fermi energies will bring infor- 
mation on the symmetry term around (below) normal density, relativistic energies 
will probe high density regions, even testing the covariant structure of the isovec- 
tor terms. The transport codes are based on mean field theories, with correlations 
included via hard nucleon-nucleon elastic and inelastic collisions and via stochastic 
forces, selfconsistently evaluated from the mean phase-space trajectory, see 13 4 5] 
Stochasticity is essential in order to get distributions as well as to allow the growth 
of dynamical instabilities. The isovector part of the EoS has been tested system- 
atically by using two different behaviors of the symmetry energy below saturation: 
one [Asysoft) where it is a smooth decreasing function towards low densities, and 
another one {Asystif f) where we have a rapid decrease, 

2. The Prompt Dipole 7- Ray Emission 

The possibility of an entrance channel bremsstrahlung dipole radiation due t o an 
initial different N/Z distribution was suggested at the beginning of the nineties 1^^^. 
After several exp erimental evidences, in fusion as well as in deep-inelastic reactions 
| ij | lU | ll | lz | lt5 | 14 | have now a good understanding of the process and stimulating 
new perspectives from the use of radioactive beams. 
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During the charge equihbration process taking place in the first stages of dissi- 
pative reactions between coUiding ions with different N/Z ratios, a large amplitude 
dipole collective motion develops in the composite dinuclear system, the so-called 
dynamical dipole mode. This collective dipole gives rise to a prompt 7-ray emission 
which depends: i) on the absolute value of the intial dipole moment 
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being Rz = ^'^^^^ and Rn = the center of mass of protons and of neu- 

trons respectively, while Rp and Rt are the projectile and target radii; ii) on the 
fusion/deep- inelastic dynamics; iii) on the symmetry term, below saturation, that 
is acting as a restoring force. 

A detailed description is obtained in a microscopic approach based on semi- 
classical transport equations, of Landau- Vlasov type, ^1^, where me an fiel d and 
two-body collisions are treated in a self-consistent way, see details in [EIIl!. Real- 
istic effective interactions of Skyrme type are used. The resulting physical picture 
is in good agreement with quantum Time-Dependent-Hartree-Fock calculation 
In particular we can stu dy in detail how a collective dipole oscillation develops in 
the entrance channel ^1^. 

We can follow the time evolution of the dipole moment in the r-space, D{t) = 
^^{Rz — Rn) and in p— space, DK{t) — (-^ — ■^), with Pp center of mass 
in momentum space for protons (neutrons), is just the canonically conjugatemo- 
mentum of the D{t) coordinate, i.e. as operators [D{t), DK{t)\ = ih see 
A nice "spiral-correlation" clearly denotes the collective nature of the mode. 

From the dipole evolution given from the Landau - Vla sov transport we can di- 
;tly apply a bremsstrahlung ("brems 
photon emission probability {E^ = Hlu): 



rectly apply a bremsstrahlung ("bremss") approach to estimate the "direct" 



where D"{uj) is the Fourier transform of the dipole acceleration D"(t). We remark 
that in this way it is possible to evaluate, in absolute values, the corresponding 
pre-equilibrium photon emission. 

We must add a couple of comments of interest for the experimental selection 
of the Dynamical Dipole: i) The centroid is always shifted to lower energies (large 
deformation of the dinucleus); ii) A clear angular anisotropy should be present 
since the prompt mode has a definite axis of oscillation (on the reaction plane) at 
variance with the statistical GDR. In a very recent experiment the prompt dipole 
radiation has been investigated with a 47r gamma detector. A strong dipole-like 
photon angular distribution (9^) = Wq [1 + 02^2(2036'^)], 6^ being the angle between 
the emitted photon and the beam axis, has been observed, with the 02 parameter 
close to —1, seetI31. 
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Fig. 1. Upper Curves: ^^-^Sn +^®Ni system (E = WAMeV, b = 4fm). Lower Curves: same 
reactions but induced by ^^*Sn. Left Panel: DK-D Spirals. Central panel: neutron (upper) and 
proton (lower) emissions. Right panel: 7 spectrum. Asy-soft (solid lines) and asy-stiff (dashed 
lines) symmetry energies. 



The deviation from a pure dipole form can be interpreted as due to the rota- 
tion of the dinucleus symmetry axis vs. the beam axis during the Prompt Dipole 
Emission. From accurate angular distribution mesuremcnts we can then expect to 
get a direct information on the Dynamical Dipole Life Time. 

At higher beam energies we expect a decrease of the direct dipole radiation for 
two main reasons both due to the increasing importance of hard NN collisions: i) 
a larger fast nucleon emission that will equilibrate the isospin before the collective 
dipole starts up; ii) a larger damping of the collective mode due to np collisions. 

The prompt dipole radiation also represents a nice cooling mechanism on the 
fusion path. It could be a way to pass from a warm to a cold fusion in th e sy nthesis 
of heavy elements with a noticeable increase of the survival probability, 

Psurv^dipole P^Psurvi,-^ ^7) 1 /-1 r> \ 1 -^surv{-^ ^7) -1 /q\ 

^ surv ^ surv J ^ surv V-'^ / 



Symmetry Energy Effects 

The use of unstable neutron rich projectiles would largely increase the effect, due to 
the possibility of larger entrance c hann el asymmetries. In order to suggest proposals 
for the new RIB facility Spiral 2,1^1] we have studied fusion events in the reaction 
^^'^Sn+^^ Ni at lOAMeV,^. We espect a Monster Dynamical Dipole, the initial 
dipole moment D{t — 0) being of the order of 50fm, about two times the largest 
values probed so far, allowing a detailed study of the symmetry potential, below 
saturation, responsible of the restoring force of the dipole oscillation and even of 
the damping, via the fast neutron emission. In the Fig.l (upper) we present some 
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preliminary very promising results. The larger value of the symmetry energy for the 
Asysoft choice at low densities, where the prompt dipole oscillation takes place, 
leads to some clear observable effects: i) Larger Yields, as we see from the larger 
amplitude of the "Spiral" (left panel) and finally in the spectra (right panel); ii) 
Larger mean gamma energies, shift of the centroid to higher values in the spectral 
distribution (right panel); iii) Larger width of the "resonance" (right panel) due to 
the larger fast neutron emission (central panel). We note the opposite effect of the 
Asy-stiffness on neutron vs proton emissions. The latter point is important even for 
the possibility of an independent test just measuring the N/Z of the pre-equilibrium 
nucleon emission. 

The symmetry energy influence can be of the order of 20%, and so well detected. 
In the lower part of the same figure we present the same results for reactions induced 
by the stable ^'^'^'Sn beam. We still se the Iso — EoS effects, but largely reduced. 



3. Isospin Dynamics in Neck Fragmentation at Fermi Energies 

It is now quite well established that the largest part of the reaction cross section for 
dissipative collisions at Fermi energies goes through the Neck Fragmentation chan- 
nel, with IMFs directly prod uced in the interacting zone in semiperipheral colli- 
sions on very short time scales 12 31^4 1 -y^g predict interesting isospin transport 
effects for this new fragmentation mechanism since clusters are formed still in a di- 
lute asymmetric matter but always in contact with the regions of the projectile-like 
and target-like remnants almost at normal densities. Since the difference between 
local neutron-proton chemical potentials is given by /j.„ — /^p = ^Esym{p?,l p), we ex- 
pect a larg er neu tron flow to the neck clusters for a stiffer symmetry energy around 
saturation, [II2S1_ xhe isospin dynamics can be directly extracted from correlations 
between N / Z , alignement and emission times of the IMFs. The alignment between 
PLF — IMF and PLF — TLF directions represents a very convincing evidence of 
the dynamical origin of the mid-rapidity fragments produced on short time scales 
The form of the ^piane distributions (centroid and width) can give a direct 
information on the fragmentation mechanism ' . Recent calculations confirm that 

the light fragments are emitted first, a general feature expected for that rupture 
oh 

mechanism . The same conclusion can be derived from direct emission time mea- 
surements based on deviations from Viola systematics observed in event- by-event 
velocity correlations between IMFs and the PLF /TLF residues I^SElESI ^^^^ 
figure out a continuous transition from fast produced fragments via neck instabili- 
ties to clusters formed in a dynamical fission of the projectile(target) residues up to 
the evaporated ones (statistical fission) . Along this line it would be even possible to 
disentangle the effects of volume and shape instabilities. A neutron enrichment of 
the overlap ("neck") region is expected, due to the neutron migration from higher 
(spectator) to lower (neck) density regions, directly related to the slope of the sym- 
metry energy l^Sl, A very nice new analysis has been presented on the Sn + Ni data 
at 35 AMeV by the Chimera CoUab., Fig. 2 of ref.l^. A strong correlation between 
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neutron enrichemnt and alignement (when the short emission time selection is en- 
forced) is seen, that can be reproduced only with a stiff behavior of the symmetry 
energy. This is the first clear evidence in favor of a relatively large slope ( symmetry 
pressure) around saturation. 

4. Isospin Distillation with Radial Flow 

In central collisions at 30-50 MeV/A, where the full disassembly of the system 
into many fragments is observed, one can study speci fically pro perties of liquid- 
gas phase transitions occurring in asymmetric matter EEHSH] Por instance, in 
neutron-rich matter, phase co-existence leads to a different asymmetry in the liquid 
and gaseous phase: fragments (liquid) appear more symmetric with respect to the 
initial matter, while light particles (gas) are more neutron-rich. The amplitude 
of this effect depends on specific properties of the isovector part of the nuclear 
interaction, namely on the value and the derivative of the symmetry energy at low 
density. 

This investigation is interesting in a more general context: In heavy ion collisions 
the dilute phase appears during the expansion of the interacting matter. Thus we 
study effects of the coupling of expansion, fragmentation and distillation in a two- 
component (neutron-proton) system. We present fully consistent dynamical results, 
based on stochastic microscopic transport approaches widely tested in heavy ion 
collisions ISUISM'SS' r^^ie correct treatment of fluctuations is essential to reproduce 
the dynamics of spinodal instabilities in the dilute expansion phase. 

We focus on central collisions, b = 2 fm, considering symmetric reactions be- 
tween systems having three different initial asymmetry: ^^"^Sn -i-H^ gfi^'^i _|_124 
571,132 5-71 -^132 5-^^ with {N/Z),„ = 1.24,1.48,1.64, respectively. The considered 
beam energy is 50 MeV/A. 1200 events have been run for each reaction and for 
each of the two symmetry energies adopted (asy-soft and asystiff, see before) 

In central collisions, after the initial coUisional shock, the system expands and 
breaks up into many pieces, due to the development of volume (spinodal) and 
surface instabilities. The formation of a bubble-like configuration is observed, where 
the initial fragments are located. 

First, let us briefly recall some general features concerning the isotopic content 
of fragments and emitted nucleons, as obtained with the two considered iso-EOS's. 
In the following we will restrict our analysis to fragments with charge in the range 
between 3 and 10 (that we call intermediate mass fragments (IMF's) ). The average 
N/Z of emitted nucleons (gas phase) and IMF's is presented in Figl2]as a function of 
the initial {N/Z)i„ of the three colliding systems. One observes that, generally, the 
gas phase (right panel) is more neutron-rich in the asysoft case, while IMF's (left 
panel) are more symmetric. This is due to the larger value of the symmetry energy 
at low density for the asysoft parameterization'^'. It is interesting to note that, 
in the asystiff case, due to the rather low value of the symmetry energy. Coulomb 
effects dominate and the N/Z of the gas phase becomes lower than that for IMF's, 
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Fig. 2. The N/Z of the liquid (left) and of the gas (right) phase is displayed as a function of the 
system initial N/Z. Full lines and symbols refer to the asystiff parameterization. Dashed lines and 
open symbols are for asysoft. 



because protons are preferentially emitted. Now we move to investigate in more 
detail correlations between fragment isotopic content and kinematical properties. 
The idea in this investigation is that fragmentation originates from the break-up 
of a composite source that expands with a given velocity field. Since neutrons and 
protons experience different forces, one may expect a different radial flow for the 
two species. In this case, the N/Z composition of the source would not be uniform, 
but would depend on the radial distance from the center or mass or, equivalently, on 
the local velocity. This trend should then be reflected in a clear correlation between 
isospin content and kinetic energy of the formed IMF's, 



5. Isospin Transport in Peripheral Collisions at the Fermi Energies 

Now we investigate peripheral collisions between similar systems with different 
isospin, specifically collisions of different combinations of Sn isotopes. The isospin 
is used here both as a tracer of the reaction mechanism, as well as an observable 
of interest with respect to the iso-EOS. In the prevoius sections we have seen the 
isospin dynamics in instability regions and fragment formation. Here we study more 
directly the isospin transport in binary events. A neck of density below normal den- 
sity develops between the two heavy residues, the evolution of which is driven by 
the motion of the spectators. During this phase isospin is transferred to the neck 
due to the density difference between the neck and the residues; this effect is called 
isospin migration, which leads to a more neutron-rich neck. In addition in collision 
systems with different asymmetry isospin is transported through neck due to the 
asymmetry gradient leading to an equilibration of the isospin of the residues, which 
has been called isospin diffusion. Thus in asymmetric systems there is a competion 
of isospin migration and diffusion. 

The isospin transport can be discussed in a compact way by means of the chem- 
ical potentials for protons and neutrons as a function of density p and isospin I^^. 
From this the p/n currents can be expressed as 

jp/n=^p/„Vp-i?^/„V/ (4) 

with D'' , the drift, and D^, the diffusion coefhcients for transport, which are 
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given explicitely in ref. Of interest here are the differences of currents between 
protons and neutrons which have a simple relation to the density dependence of 
the symmetry energy 

DP - DP 4/ ^^"^" 

Di-D^^x 4pE,yrn . (5) 

Thus the isospin transport due to density gradients, i.e. isospin migration, depends 
on the slope of the symmetry energy, or the symmetry pressure, while the transport 
due to isospin concentration gradients, i.e. isospin diffusion, depends on the absolute 
value of the symmetry energy. In peripheral collisions discussed here, residues of 
about normal density are in contact with the neck region of density below saturation 
density. In this region of density a stiff iso-EOS has a smaller value but a larger 
slope compared to a soft iso-EOS. Correspondingly we expect opposite effects of 
these models on the migration and diffusion of isospin. 

We will in particular discuss the so-called Imbalance Ratio (also called Rami- 
or transport-ratio), which is defined as^^ 

^^'^ ~ (xH - x^) ■ 

Here, x is an isospin sensitive quantity that is to be investigated with respect to 
equilibration, as e.g. the asymmetry / = [N — Z)/ {N + Z) which is considered 
here primarily, but also other quantities, such as isoscaling coefficients, ratios of 
production of light fragments, etc. The indices H and L refer to the symmetric 
reaction between the heavy (n-rich) and the light (n-poor) systems, while M refers 
to the mixed reaction. P, T denote the rapidity region, in which this quantity is 
measured, in particular the projectile and target rapidity regions. Clearly, this ratio 
is ±1 in the projectile and target region, respectively, for complete transparency, 
and oppositely for complete rebound, while it is zero for complete equilibration. 



Correlation with kinetic energy loss 

The centrality dependence of the Imbalance Ratio, fo r (Sn, Sn) collisions, has been 
investigated in experiments as well as in theory l ^5 | 38 | -^ji^jj j-ather promising 
results for the sensitivity to the symmetry term stiffness. We propose here a new 
analysis which appears experimentally more selective. The i nterac tion time cer- 
tainly influences the amount of isospin equilibration, see refs. I ^^ l ^^ l. On the other 
hand, longer interaction times should be correlated to a larger dissipation. The dis- 
sipation, in turn, has been measured, i.e. in deep inelastic collisions, by the kinetic 
energy loss. 

In Fig. 3 we present the imbalance ratio as a function of the relative energy 
loss (H:124Sn, L:112Sn, two beam energies 35 and 50AMeV). In the upper panel 
we separate the results for Momentum Dependent and Independent (MD vs MI) 



I 



I 
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Fig. 3. Imbalance ratio as a function of relative energy loss: (upper panel) Separately for stiff 
(solid) and soft (dashed) iso-EOS, and for MD (circles and diamonds) and MI (squares and 
tringles)interactions, in the projectile region (full symbols) and the target region (open symbols), 
(lower panel) Quadratic fit to all points for the stiff (solid), resp. soft (dashed) iso-EOS. 

interactions and stiff and soft iso-EOS, in each case the results for 35 and 50 MeV 
are collected together and connected by lines. We now see that the points for the 
different incident energies approximately fall on one line, however with some scatter 
around it. 

The curves for the soft EOS (dashed) are generally lower in the projectile region 
(and oppositely for the target region), i.e. show more equilibration, that those for 
the stiff EOS. In order to emphasize this trend we have, in the lower panel of the 
figure, collected together the values for the stiff (circles, solid) and the soft (squares, 
dashed) iso-EOS, and fitted them by a quadratic curve. It is seen that this fit gives 
a good representation of the trend of the results. 

The difference between the curves for the stiff and soft iso-EOS in the lower panel 
then isolates the effect of the iso-EOS from kinematical effects depending on the 
interaction time. It is seen, that there is a systematic effect of the symmetry energy 
of the order of about 20 percent, which should be measurable. The correlation 
suggested in Fig. 3 should represent a general feature of isospin diffusion, and it 
would be of great interest to verify it experimentally. 

6. Relativistic Collisions 

Finally we focus our attention on relativistic heavy ion collisions, that provide 
a unique terrestrial opportunity to probe the in-medium nuclear interaction at 
high densities. An effective Lagrangian approach to the hadron interacting system 
is extended to the isospin degree of freedom: within the same frame equilibrium 
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properties {EoS, and transport dynamics 1^11421 ^.g^j-^ ]-,g consistently derived. 
Within a covariant picture of the nuclear mean field, for the description of the 
symmetry energy at saturation (a) only the Lorentz vector p mesonic field, a nd (b) 
both, the vector p (repulsive) and scalar S (attractive) effective fields I43|44 | 
included. In the latter case the competi tion b etween scalar and vector fields leads to 
a stifFer symmetry term at high density The presence of the hadronic medium 
leads to effective masses and momenta M* — M + Es, k*'^ = fc'^ — E'^, with S^, S'^ 
scalar and vector self-energies. For asymmetric matter the self-energies are different 
for protons and neutrons, depending on the isovector meson contributions. We will 
call the corresponding models as NLp and NLpS, respectively, and just NL the 
case without isovector interactions. 

For the more general NLp5 case the self-energies of protons and neutrons read: 

n) = ~J„a{ps) ± JsPsi, n) = f^f T /pis ' ■ ^PP^^ signs), (7) 

where ps = Psp + Psn, j" = jp + j", Ps3 = Psp ~ Psn, is = jp - jn are the total 
and isospin scalar densities and currents and fa,uj,p,s are the coupling constants of 
the variou s mesonic fields. (j{ps) is the solution of the non linear equation for the tr 
field '"--IIJ por the description of heavy ion collisions we solve the covariant transport 
equation of the Boltzmann type^^i^^ within the Relativistic Landau Vlasov {RLV) 
method, using phase-space Gaussian test particles SHl^ and applying a Monte-Carlo 
procedure for the hard hadron collisions. The collision term includes elastic and 
inelastic processes involving the production/absorption of the A(1232Mey) and 
N*{lAAOMeV) resonances as well as their decays into pion channels, SSSZl, A larger 
repulsive vector contribution to the neutron energies is given by the /O-cou pling. 
This is rapidly increasing with density when the 6 field is included SSH] 
consequence we expect a good sensitivity to the covariant structure of the isovector 
fields in nucleon emission and particle production data. Mor eover the presence of 
a Lorentz magnetic term in the relativistic transport equation I 41 | 42 | l | .^jj^j enhance 
the dynamical effects of vector fields SSI. 

Isospin Flows 

(n,p) flow differences will be directly affected. In ref transverse and elliptic flows 
results are shown for the ^^"^ Sn-\-^'^'^ Sn reaction at 1.5 AGeV {b — 6/m). The effect 
of the different structure of the isovector channel is clear. Particularly evident is the 
splitting in the highpt region of the elliptic flow. In the {p + S) dynamics the high-pt 
neutrons show a much larger squeeze — out. This is fully consistent with an early 
emission (more spectator shadowing) due to the larger p-field in the compression 
stage. 

Isospin effects on sub-threshold kaon production at intermediate energies 

Kaon production has been proven to be a reliable observable for the high density 
EoS in the isoscalar sector 1S1151J Here we show that the K^''^ production (in par- 



2, 2008 9:59 WSPC/INSTRUCTION FILE ditoro'beijingO? 




Fig. 4. Time evolution of the resonances and pions tt^'^ (left), and kaons {K~^'^ 

(right) for a central (fa = fm impact parameter) Au+Au collision at 1 AGeV incident 
energy. Transport calculation using the NL, NLp, NLpS and DDF models for the iso- 
vector part of the nuclear EoS are shown. The inset contains the differential / K'^ 
ratio as a function of the kaon emission time. 

ticul ar the K '^/K'^ yield ratio) can be also used to probe the isovector part of the 

Using our RMF transport approach we have analyzed pion and kaon produc- 
tion in central ^^'^Au +197 collisions in the 0.8 — 1.8 AGeV beam energy range, 
comparing models giving the same "soft" EoS for symmetric matter and with dif- 
ferent effective field choices for Egym- We will use three Lagrangians with constant 
nucleon- meson coupl ings {NL... type, see before) and one with density dependent 
couplings (DDF , see'^^, recently suggested for better nucleonic properties of neu- 
tron stars ESlMl, 

Fig. m reports the temporal evolution of A^'°'++ resonances, pions (tt^-") and 
kaons (_fi'+'°) for central Au-|-Au collisions at lAGeV. It is clear that, while the 
pion yield freezes out at times of the order of 50/to/c, i.e. at the final stage of 
the reaction (and at low densities), kaon production occur within the very early 
(compression) stage, and the yield saturates at around 20/m/c. Consistently, as 
shown in the insert, larger isospin effects are expected for emitted kaons. 

When isovector fields are included the symmetry potential energy in neutron- 
rich matter is repulsive for neutrons and attractive for protons. In a HIC this leads 
to a fast, pre-equilibrium, emission of neutrons. Such a mean field mechanism, 
often referred to as isospin fractionation 49il]^ jg responsible for a reduction of the 
neutron to proton ratio during the high density phase, with direct consequences on 
particle production in inelastic NN collisions. 

Threshold effects represent a more subtle point. The energy conservation in a 
hadron collision in general has to be formulated in terms of the canonical momenta, 
i.e. for a reaction 1+2 3+4 as Si„ = i^i^^'i)'^ = (^3 +''^4 )^ ~ ^out- Since hadrons 
are propagating with effective (kinetic) momenta and masses, an equivalent relation 
should be formulated starting from the effective in- medium quantities fc*'' = fc^ — E'' 
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and TO* = TO + Es, where and are the scalar and vector self-energies, Eqs.(l7l). 
The self-energy contributions will influence the particle production at the level of 
thresholds as well as of the phase space available in the final channel. 

In the few AGeV region the effect is larger for the /K^ compared to the 
7r^/7r+ ratio. This is due to the subthreshold production and to the fact that the 
isospin effect enters twice in the two-step production of kaons, see S3. Interestingly 
the Ibo-EoS effect for pions is increasing at lower energies, when approaching the 
production threshold. 

We have to note that in a study of kaon production in excited nuclear matter the 
dependence of the K° /K~^ yield ratio on the effective isovector interaction appears 
much larger (see Fig. 8 of ref.^Sl). The point is that in the non-equilibrium case of a 
heavy ion collision the asymmetry of the source where kaons are produced is in fact 
reduced by the n p "transformation" , due to the favored nn pA^ processes. 
This effe ct is almost absent at equilibrium due to the inverse transitions, see Fig. 3 
of ref.SSl. Moreover in infinite nuclear matter even the fast neutron emission is not 
present. This result clearly shows that chemical equilibrium models can lead to 
uncorrect results when used for transient states of an open system. 



7. Testing Deconfinement at High Isospin Density 

The hadronic matter is expected to undergo a phase transition into a deconfined 
phase of quarks and gluons at large densities and/or high temperatures. On very 
general grounds, the transition's critical densities are expected to depend on the 
isospin of the system, but no experimental tests of this dependence have been 
performed so far. 

Here we suggest the possibility of an earlier transition density to a mixed phase 
in isospin asymmetric systems'^. Concerning the hadronic phase, we use the rela- 
tivistic non-linear model of Glendenning-Moszkowski (in particular the "soft" GM3 
choice) where the isovector part is treated just with p meson coupling, and the 
iso-stiffer NLp6 interaction. For the quark phase we consider the MIT bag model 
with various bag pressure constants. In particular we ar e int erested in those param- 
eter sets which would allow the existence of quark stars i.e. p arameters sets for 
which the so-called Witten-Bodmer hypothesis is satisfied I^HEII One of the aim of 
our work it to show that if quark stars are indeed possible, it is then very likely to 
find signals of the formation of a mixed quark-hadron phase in intermediate-energy 
heavy-ion experiments 

Th e structure of the mixed phase is obtained by imposing the Gibbs conditions 
| 60 | 61| ciigjnica,! potentials and pressure and by requiring the conservation of the 
total baryon and isospin densities 

{H) _ (Q) {H) _ (Q) 

p(«)(T,MS)-^^^^mMS). 
PS = (1 - x)Pb + XPb ' 
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Fig. 5. Variation of the transition density with proton fraction for various hadronic EoS 
parameterizations. Dotted line: GM3 parametrization; dashed line: NLp parametrization; 
solid line: NLpS parametrization. For the quark EoS, the MIT bag model with B^^^=150 
MeV. The points represent the path followed in the interaction zone during a semi-central 
^^^Sn+l^^Sn collision at 1 AGeV (circles) and at 300 AMeV (crosses). 

P3-(l-x)pf + XP?, (8) 

where X is the fraction of quark matter in the mixed phase. In this way we get the 
binodal surface which gives the phase coexistence region in the (T, ps, pa) space. 
For a fixed value of the conserved charge ps we will study the boundaries of the 
mixed phase region in the (T, pb) plane. In the hadronic phase the charge chemical 
potential is given by /is ~ 2Esym{pB)j^ ■ Thus, we expect critical densities rather 
sensitive to the isovector channel in the hadronic EoS. 

In Fig. [5] we show the crossing density pcr separating nuclear matter from the 
quark-nucleon mixed phase, as a function of the proton fraction Z/A. We can see 
the effect of the (5-coupling towards an earlier crossing due to the larger symme- 
try repulsion at high baryon densities. In the same figure we report the paths in 
the (p, Z/A) plane followed in the cm. region during the collision of the n-rich 
"'^^^Sn-^-'^^^Sn system, at different energies. At 300 AMeV we are just reaching the 
border of the mixed phase, and we are well inside it at 1 AGeV . We expect a neu- 
tron trapping effect, supported by statistical fluctuations as well as by a symmetry 
energy difference in the two phases. In fact while in the hadron phase we have a 
large neutron potential repulsion (in particular in the NLpS case), in the quark 
phase we only have the much smaller kinetic contribution. Observables related to 
such neutron "trapping" could be an inversion in the trend of the formation of neu- 
tron rich fragments and/or of the 7r^/7r+, / yield ratios for reaction products 
coming from high density regions, i.e. with large transverse momenta. In general we 
would expect a modification of the rapidity distribution of the emitted "isospin" , 
with an enhancement at mid-rapidity joint to large event by event fluctuations.. 

8. Perspectives 

We have shown that violent collisions of n-rich heavy ions from low to relativistic 
energies can bring new information on the isovector part of the in-medium inter- 
action, qualitatively different from equilibrium EoS properties. We have presented 



2, 2008 9:59 WSPC/INSTRUCTION FILE ditoro'beijingO? 



quantitative results in a wide range of beam energies. At low energies we see isospin 
effects on the dissipation in deep inelastic collisions, at Fermi energies the Iso-EoS 
sensitivity of the isospin transport in fragment reactions and finally at intermediate 
the dependence of differential flows on the Iso — MD and effective mass splitting. 
In relativistic collisions we have shown the possibility of a direct measure of the 
Lorentz structure of the isovector fields at high baryon density. We have presented 
quantitative results for differential collective flows and yields of charged pion and 
kaon ratios. Important non-equilibrium effects for particle production are stressed. 
Finally our study supports the possibility of observing precursor signals of the phase 
transition to a mixed hadron-quark matter at high baryon density in the collision, 
central or semi-central, of neutron-rich heavy ions in the energy range of a few GeV 
per nucleon. In conclusion the results presented in this "Isospin Journey" appear 
very promising for the possibility of exciting new results from dissipative collisions 
with radioactive beams. 

Acknowledgements. We warmly thanks A.Drago and A.Lavagno for the intense col- 
laboration on the mixed hadron-quark phase transition at high baryon and isospin 
density. 
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